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Abstract 

As a step toward clarification of tlie power of supersymmetry (SUSY) in Matrix the- 
ory, a complete calculation, including all the spin effects, is performed of the effective 
action of a probe D-particle, moving along an arbitrary trajectory in interaction with a 
large number of coincident source D-particles, at one loop at order 4 in the derivative 
expansion. Furthermore, exploiting the SUSY Ward identity developed previously, the 
quantum-corrected effective supersymmetry transformation laws are obtained explicitly 
to the relevant order and are used to verify the SUSY-invariance of the effective action. 
Assuming that the agreement with 11-dimensional supergravity persists, our result can 
be regarded as a prediction for supergravity calculation, which, yet unavailable, is known 
to be highly non-trivial. 

'l'kazaina@hep3.c.u-tokyo.ac.jp 
^ t etsu@hep 1 . c . u-tokyo .ac.jp 



1 Introduction 



One of the main themes of string theory in the past few years is the correspondence 
between supergravity/ superstring theories in the bulk and certain types of supersymmetric 
theories defined on the worldvolume of various brane configurations or "on the boundary" 
of the space-time. The proposal of Maldacena|jl|, termed AdS/CFT correspondence, has 
by now been extended to much wider class of systems than was originally envisaged and 
has spurred wide variety of new developments 0. 

Although conjectured prior to this proposal, the Matrix theory for M-theory, put for- 
ward by Banks, Fischler, Shenker and Susskind , can be thought of as a prime example 
of this correspondence. Re- interpreted in the frame work of discrete light-cone quantiza- 
tion [0, it has enjoyed numerous successes. Just to mention only the direct comparison 
with eleven dimensional supergravity, complete agreement for the multi-graviton scat- 
tering (including the recoil effects) at 2-loop0, ^ and that for the two-body potential 
between arbitrary fermionic as well as bosonic objects at 1-loop should be regarded as 
highly non-trivial and remarkable. 

Just as the mechanism of the general bulk-boundary correspondence has not been 
clearly identified, the origin of these successes in the Matrix theory is yet to be fully 
understood. Now in string theory, it is often the case that symmetry principles play 
decisive roles, much more so than in local field theories, in determining the dynamics of 
the system, and the situation should be the same in Matrix theory, which is an explicit 
representation of M-theory that unifies all the string theories. 

Indeed there have been a number of studies p, ^, |Tl|, |12|, |1^ which point to the 
assertion that, in particular, the high degree of supersymmetry, namely with the maxi- 
mally allowed 16 supercharges, is powerful enough to determine the effective action of the 
D-particles, which can be directly compared with the corresponding supergravity calcula- 
tion. If it is indeed the case, it is rather surprising since usually such a global symmetry 
can only give certain relations among the correlation functions and not more, and con- 
sequently the dynamical significance of Matrix theory would have to be reconsidered. 
However, upon close examinations, even restricted to the simple so-called "source-probe" 
situation, one can argue that the existing analysis is not quite complete. This is largely 
due to the fact that for a system with maximal supersymmetry unconstrained superfield 
formulation does not exist and therefore a clear-cut off-shell analysis is not possible: One 
is forced to deal with the component formalism, where the supersymmetry algebra gets 
intertwined with gauge symmetry and does not close without the aid of equations of 
motion. 
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In a previous work[n|, we gave a rather extensive discussion on the existing htera- 
ture and emphasized the importance of the consistency of the approximation scheme in 
studying the symmetry of the effective action, which is an off-shell entity. We argued that 
the only consistent procedure is to deal with the trajectory, including the spin degrees of 
freedom, with arbitrary time dependence and adopt the derivative expansion according 
to the concept of order, defined as the number of derivatives plus half the number of 
fermions. Only this scheme is free of total derivative ambiguities inevitably present in the 



effective action. For more discussions, see the Section 3 of 14 



Following this philosophy, we derived in |jT^ a completely off-shell SUSY Ward identity 
in the background gauge (which is naturally intertwined with the BRST symmetry) and 
applied it to the effective action at order 2 to study the power of supersymmetry. Our 
conclusion was that at this order the effective action is indeed determined by the Ward 
identity alone. As we already remarked there, however, such a result was to be expected 
since at order 2 the higher derivatives, such as the acceleration etc., can be eliminated from 
the effective action by integration by parts and our analysis was essentially the same as in 
the existing literature. The full significance of the completely off-shell analysis becomes 
apparent starting from the next order, i.e. from order 4, where complete elimination of 
higher derivatives is no longer possible. 

Unlike the case of order 2, at which the 1-loop effective action was already available, the 



only fully off-shell results so far known at order 4 are the famous bosonic part 1 15, 16, 17 



given (in Euclidean formulation) by —N J dr {15 / 16)v'^ / r"^ , with the number of source 
D-particles, and a part of the fermionic contributions containing 8 powers of the fermion 
field ^ [1^. The latter is the simplest among the fermionic contributions since it does not 
involve any derivatives. A calculation of 0{9'^) contribution, which is already formidable, 
was attempted in |]l9l but was not fully completed. As for 0{6^) and 0{6^) contributions, 
no attempt has been made to date. 

In this work, as a necessary step toward clarification of the power of supersymmetry 
in Matrix theory, we perform a complete calculation of the effective action for a probe D- 
particle, including all the spin effects, at one-loop at order 4, with the aid of the algebraic 
manipulation program Mathematica. Furthermore, exploiting the SUSY Ward identity 
developed previously, the quantum-corrected effective supersymmetry transformation laws 
are obtained explicitly to the relevant order and are used to verify the SUSY-invariance 
of the effective action so obtained. Since the corresponding supergravity calculation is 
available only up to 0{9'^) we cannot at present time test if our result agrees with 
supergravity. Rather, provided that the agreement with supergravity persists, our result 
should be regarded as a prediction until such a calculation will have been made. 
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Organization of the rest of this article is as follows: We start in Section 2 with a brief 
summary of the Matrix theory and its symmetries, mainly to set the notations. Section 
3 is devoted to the calculation of the effective action at one-loop. After explaining the 
derivative expansion scheme used and reviewing the existing results for the effective ac- 
tion in Section 3.1, we describe the calculational procedures in Section 3.2. The result of 
the calculation, together with how it was simplified using SO (9) Fierz identities, is pre- 
sented in Section 3.3. Since the actual calculation, performed by Mathematica, involved 
an enormous number of steps, it is desirable to have an independent check. For this 
purpose, as well as for its own interest, we compute in Section 4 the quantum-corrected 
SUSY transformation laws and check the invariance of the effective action under these 
transformations: Following a brief review of the SUSY Ward identity in Section 4.1, we 
sketch the calculational procedure in Section 4.2. Then in Section 4.3 we describe how 
the invariance under SUSY was verified. Finally, we give a summary and discussions in 
Section 5. Three appendices, A~C are provided to supply some details of the calculations. 
In appendix A we describe a new efficient algorithm for generating SO (9) Fierz identities. 
Non-trivial two-point functions needed for the calculation of the SUSY transformation 
laws are collected in Appendix B and the quantum-corrected SUSY transformation laws 
are displayed in Appendix C. 

2 Preliminaries 

We begin with a very brief summary of the Matrix theory and its symmetries, mainly to 
set our notations. 

The classical action for the U{N + 1) Matrix theory in the Euclidean formulation is 
given by 



In this expression, X^{t), Aij{T) and Oa,jj(r) are the {N -\- 1) x (A^+ 1) hermitian matrix 
fields, representing the bosonic part of the D-particles, the gauge fields, and the fermionic 
part of the D-particles, respectively. D^- = dr — igA is the covariant derivative, 7"^ are 
the real symmetric 16 x 16 5*0(9) 7-matrices, and the vector index m runs from 1 to 9. 
We put a tilde on some relevant symbols to remind us of Euclidean formulation. 

This action is known to possess a number of important symmetries. The first is the 




(2.1) 
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obvious global Spin(9) invariance, which nevertheless is quite non-trivial in the fermionic 
sector being responsible for various Fierz identities to be used extensively later. The 
second is the CPT symmetry inherited from the 10- dimensional super Yang-Mills theory 
from which the above action can be obtained by dimensional reduction. The third is the 
invariance under the U{N + 1) gauge transformations given, with the gauge parameter 
matrix A, by 

6aA = [D^, A] , 5aX^ = ig [A, X^] , 

6Ae = tg[A,&] . (2.2) 

The fourth, and the main focus of this article, is the supersymmetry with 16 spinorial 
parameters ea- The transformation laws are 

= e^e , ^.X™ = -ie^-f^'Q , (2.3) 
5,e = t {[D,, X^] 7- + I X„] 7™") e . (2.4) 

Although the algebra closes only on-shell up to field-dependent gauge transformations, Sq 
itself is invariant without the use of equations of motion, i. e. off-shell. 

In addition to these well-known symmetries, there is a so-called generalized conformal 



symmetry [pT|, p^ , p3| , which may be used to restrict the form of the effective action. 
Finally, the agreement with the 11-dimensional supergravity calculation for the multi- 
body processes]^] strongly suggests that the 11-dimensional Lorentz symmetry is hidden 
in 5*0, awaiting to be disclosed. 

In this article, we shall concentrate on the so-called source-probe situation, namely 
the configuration of a probe D-particle interacting with a large number, N, of the source 
D-particles all sitting at the origin. This is expressed by the splitting 

Xmir) = -B^r) + YUr) , e„(r) = -e^{r) + ^„(r) , (2.5) 

9 9 

N N 

Bm{r) = diag(r^(r),M7^^) , ^^(r) = diag (^^„(r), oX^~To) , (2.6) 

where BmiT) and ^a(r) are the bosonic and the fermionic backgrounds expressing the 
positions and the spin degrees of freedom of the D-particles respectively and Ymir) and 
\l/a(T) denote the quantum fluctuations around them. We will be interested in the general 
case where rm{T) and 6'a(r) are arbitrary functions ofr not satisfying equations of motion. 

In order to quantize the system and perform the calculation of the effective action of 
the probe D-particle, we must fix the gauge. Practically the only tractable choice, and 
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indeed the one used for all the calculations in the past, is the background gauge specified 
by the gauge-fixing function G given by 

G = -d^A + i[B'^,Xj^ . (2.7) 

The associated BRST transformations for the quantum part of the fields are given by 

5bA = [Dr, C] , (5bF„ = -tg C] , 

6B^ = zg{C,Q} , (2.8) 
6bC = igC^ , 5bC = ib , (5b6 = , 

where C, C and b are, respectively, the ghost, the anti-ghost and the Nakanishi-Lautrup 
auxiliary field. As usual, the gauge-fixing and the ghost part of the action are given 
altogether by the BRST total variation 



^gg = (^eTr / dr 



(2.9) 



3 Calculation of the OfF-shell Effective Action 



We are now ready to start the computation of the off-shell effective action for the probe 
D-particle. 

3.1 Derivative expansion and the existing results 



As was already emphasized in the introduction and further elaborated in , it is impor- 
tant that the approximation scheme for computing the effective action must be consistent 
with the freedom of adding total derivatives. The only such scheme is the derivative 
expansion according to the concept of order defined as 



order = number of r-derivatives + -number of fermions . (3.10) 

In other words, we assign order(r) = , order (9) = 1, and order (6') = 1/2, where d denotes 
the derivative with respect to r. 

This concept can be applied loop-wise. For instance, the tree level action for and 9^ 
is of the form (we use Vm to denote r^, the superscript in parenthesis is the loop number 
and the subscript signifies the order) 
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and is entirely of order 2. The order 2 part at 1 loop was first computed[] in 
given by 



dr 



and is 



(3.12) 



Consider now the order 4 part at 1-loop, which consists of the bosonic and the fermionic 



parts. The former was computed long ago in eikonal approximation in [15| and more 



recently for fully off-shell case in [16, 17] and has the well-known form 



(4),b 



-A^ dr . 

16 r7 



(3.13) 



As we shall see shortly, calculation of the fermionic part 0, which can be further classified 
by the (even) number of 6''s up to 8, is exceedingly more difficult. The 0{d^) part is 
relatively easy to compute |18| since it cannot contain any derivatives. The calculation 
of the (9(6'^) part (with three (9's) was attempted in |19] but was not fully completed. As 
for the 0{d^) and 0{6^) contributions, even an attempt has not been made. 

In the following, we perform the complete calculation of all the fermionic order 4 terms 
at 1 loop. This is made possible with the extensive use of Mathematica. 



3.2 Calculational procedure 

As we shall perform the calculation at 1 loop, we only need the part or the gauge-fixed 
action quadratic in the quantum fluctuations. To simplify the presentation as well as the 
subsequent calculations, define the (1 + 9)-component fermionic and bosonic vectors in 
the following way: 




(3.14) 



where the matrix indices are suppressed. Also define the kinetic operators and Dp for 
bosonic and fermionic fields respectively as 

- ( 2»» A-m1 ) • + <''^'' 

Here, A^^ = —9^ + r^(r) is the basic kinetic operator, and its inverse A, the basic 
propagator, will appear frequently in the actual calculations. Note that the "mass" r(r) 
is an arbitrary function of r and this will make the computation non-trivial. 
^Here and hereafter, all the results are in the background gauge. 

^Within the eikonal-type approximation, spin effects have been discussed in[^, HI H^, HI H^, |30| . 
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With these notations, the quadratic part of the action can be written as 

1 



5(2) = dTj2 



1 -.T",T, , l,Tr ■='T/ 



(3.16) 



Performing the Gaussian integration, we obtain the formal expression for the 1-loop ef- 
fective action T^^^ consisting of the bosonic and the fermionic parts: 

rw = rjj') + f (3.17) 

f[J^ = trln(l - 4w"Av"A) - ^Tr ln(l + ^A) , 

f « = -Trln (l - D'/E^D^'e) = ^r {D^'E^D^'eY ^^^^^^ 

n=l 

In the above, "tr" is the trace over the matrix indices and over the function space, while 
"Tr" further includes the trace over the spinor indices as well. At order 4, the bosonic 
part Fg'* reproduces the known result ( p.l3| ). Our interest is in the fermionic part Fp^ 
and we need to compute up to n = 4. 

Now to proceed, we must expand the propagators D^^ and Dp^ in powers of the 
derivatives. The explicit form of D^^ is given by 



1 _ / A(l - Av^Av^A)-^ 2zA(l - Av^Av^A^'^v'^A 
B - 1^ -2iAw"A(l - Av^Av^A)-'^ A{6nm - 4w"Aw"'A)-i 



(3.19) 



where the expansion of the following type should be substituted: 

{6nm - Av^'Av^'A)-^ = 6nm + Av^'Av'^A + 16t;" Ai;^ At;^ At;" A + ■ ■ ■ . (3.20) 

The corresponding expansion for Dp is given by 

D-i = {d + f)-' = -{d - /)(1 + A|^)-^A 

= -(9-/)A + (9-/)A|^A + --- . (3.21) 

When these expansions are implemented, evidently the expressions for Tr (Z^p^S^Dg^S)" 
get rather involved. Below we only display the intermediate result (which still has to be 
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fully expanded) for the simplest case of n = 1 as an example: 



Tr 



X 



- (a - /)(1 - + A|^A|^ - A|^A|^A|^)„,A 

- e.Ae^ - Ae.Av^Av^Ae^ - 2Ae,^sa^0s - SAe.v^Av^A^sa^Os 
+ 20^A|^^,A^, + 89^Af,^,Av'Av'Ae^ + O^Y^^Aes-f^^ 

+ 49^Af,.,,A^sM)] + 0{d^) . (3.22) 



The number of terms are much larger for n > 2 and altogether they add up to about 
1000. 

The next step is the actual evaluation of the trace "Tr" , which here means tracing over 
the function space and over the spinor indices. While the latter, which is nothing but the 
7-matrix algebra, is conceptually simple, the former step requires some explanation. To 
perform it efficiently, we employ the so-called "normal-ordering method", an algorithm 
invented by Okawa in The idea is to bring the expressions involving /(r) (some 



function of r), d and A in various orders into a "normal-ordered" form /(r)9™A" by 
recursively using the formulas 

df = fd + f, 

A/ = /A + A(/ + 2/9)A, (3.23) 
A9 = 9A + 2A(r ■r)A. 

Further, by using the trivial relation d'^A = r^A — 1, which follows from the definition of 
A, one can reduce the number of 9's in the normal-ordered form down to either zero or one. 
Once each term is brought to this standard form, the remaining task is to compute the 
matrix elements (r|A"|r') and (r|9A"|r'), where the latter is actually expressed in terms 
of the former as (r|(9A"|r') = |9(r| A^jr'). Now employing an integral representation, 
one can write 

(r|A"|r') = -—^——-S{t - t') 

= , / da cr"-ie-'^(-^'+^^(")')5(r - r') . (3.24) 

\n - 1)! io 

Since the exponent of the exponential in the integrand contains the operator d'^ and the 
function r^(r) which are non-commuting, one needs to make use of the Baker-Campbell- 
Hausdorff formula up to the appropriate order in the derivative expansion, so that e°"^^ 



acts directly on 5(r — r'). For more details, see |1T7| , An example of the result of such 
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a calculation is 

^rl^^^r)^-^.'-}^-l^-1^.0ia^). (3.25) 

Sir' Zobr^-^ 128 r^^ 1/8 r^^ 

3.3 Simplification procedure and the result 

As already mentioned, the procedures described above were executed by developing an 
elaborate codes for Mathematica. The number of terms at the input stage is about 1000. 
As one performs the normal-ordering and the 7-matrix algebra, this number remains to 
be roughly of the same order. We then tried to simplify the results by bringing them 
to appropriate standard forms via integration by parts and identification of the same 
structures with different repeated indices. This manipulation brought the number down 
to about 70. 

The final step is to simplify them as much as possible via the use of 5*0(9) Fierz 
identities. According to the terminology of 0], a Fierz identity in which n of the indices 
of the 7- matrices involved are not contracted among themselves is called an "n-free- index" 
identity. It turns out that, including the situations that occur in the next section where we 
examine the SUSY invariance of the effective action, we need up to 5-free-index identities 
in this parlance. Although some class of Fierz identities are known in the literature and 
a general procedure to generate all possible identities was described in |^, this was not 
enough for our purposes. The reason is two-fold: First, since we deal with 6{t) with 
arbitrary r dependence, 6, 6, 6, etc. must be treated as different spinors and hence we 
need general forms of the identities. The ones available in the literature do not cover such 
cases in sufficient generality. Second, the algorithm of turned out to be prohibitively 
time-consuming when the number of free indices exceeds 3. To overcome this difficulty, 
we developed a new more efficient algorithm, which is described in the Appendix A. None 
the less, one must cope with the fact that the number of independent Fierz identities 
grows like ~ 5 x 2" and besides the length of each identity increases rapidly with n as 
well. Consequently, it was a difficult task to find the right identities to be applied for 
simplification. 

However, when the right identities were applied, the result became remarkably simple. 
The following, we believe, are the simplest forms for the desired fermionic part of the 
1-loop effective action at order 4: 
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rw = dr 



25 (69) 35 {r-vf (99) 5 (r ■ a) 



16r7 8r9 Ar'^ 



, 105v^r,v,{9Y^9) , 5 1;, a, (^7'^'^) ^^ 

+ 4^5 32^9 + 16^7 J' ^^-^^^ 

f (1) - [,i ( 25 (^^) (^^) 245 Ti vj {99) {9Y^9) 105 Vi vj {9^'^9) {9-f^''9) 

35 u aj {9Y''9) {9-i^^9) 945 Vk vi {9Y''9) {9-f^^9) 5 (9^9) (9^9) 
12879 128 ^ 327^ 

35 n Vj {9-1^^9) {9Y''9) 35 vj {9Y^9) {9^^^9) \ 
+ ^479 ^^79 J' (3.27) 

f(i) _ [a f 21 urj (99) {9Y''9){9-fi''9) 105 riVj {9 9) {9-f^'' 9) {9-f''^9) 

1155 ri Tj Tk vi {9Y"'9) {9-f^"'9) {9j''^9) 21 rj {9-f^^9) {9-f''^9) {9Y^9) 
256713 64 

(3.28 

_ fr f 15 {9Y^9) {9-i^^9) {9-i^^9) {9-i^^9) 165 r^r^- {9-i^^9) {9-f^'^9) {9-i^^9) {9y"'9) 
1 os — / dT 



1024 512 ri3 

2145 ri Tj Tfc ri {9Y"'9) {9-f^"'9) {9-f^''9) {9y''9) ' 



2048 r 



15 



(3.29) 



Here and for the rest of the article, we omit the overall common factor of for simplicity. 
This constitutes one of the main results of this work. We remark that F^a'' is slightly 
different from the one quoted in a previous attempt while , which does not contain 



any derivatives, agrees completely with the calculation of ||T^. Furthermore, if we drop 
the terms containing the derivatives of 9 and the acceleration a, our result coincides (up 



to an overall constant) with the previous calculations |2^, |2^, |TT| for this very special 
configuration. As for comparison with 11-dimensional supergravity, we shall make a 
remark in the final section. 

In view of the fact that the calculation consisted of enormous number of steps and 
that we deal with fully off-shell configurations, the agreement cited above for special cases 
should only be regarded as a strong but not decisive evidence for the correctness of our 
calculation: It is certainly desirable to have an independent check. In the next section, 
we shall perform such a test by computing the effective SUSY transformation laws for the 
effective action and demonstrate that indeed the effective action shown above is invariant 
under these transformations. 
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4 Effective SUSY Transformation Laws and Invari- 
ance of the Effective Action 



4.1 A brief review of the SUSY Ward identity 

Since the calculation of the effective SUSY transformation laws will be based on the 
super symmetric Ward identity in the background gauge derived in a previous publication 
jl^ , let us briefly review its salient features. 

The derivation of the SUSY Ward identity in the path-integral formalism is more or 
less a textbook matter, except for two non-trivialities. One stems from the intertwining of 
the SUSY and BRST symmetries. In order to bring the SUSY Ward identity in a useful 
form, one must make judicious uses of the BRST Ward identities. The other complication 
is due to the two different origins of the Bm dependence, one from the separation of the 
background and the quantum fluctuation, namely Xm = {l/g)Bm + Ym, and the other 
from the gauge-fixing function G = —drA + i [B"^,Xm\- In order to obtain the correct 
Ward identity, one must carefully distinguish these two origins. For more details, see 

0- 

Through the procedure outlined above, one obtains the SUSY Ward identity in the 
form where the effective SUSY transformation laws can be read off in closed forms. 
Adapted to the source-probe situation under consideration, it takes the form 

where the effective SUSY transformation laws are given by 

Aer„(r) = jdr' T^%{r' ,r){{6Sn{T')) + {6Byn{r'W) , (4.31) 

- jdr'dr" T^U^",r'){6BMr)On{r')) 

X {SSUr") + SBym{r")0,) ■ (4.32) 

In the expressions above, 6^ is the SUSY variation (|2.3|), 6b is the BRST variation 



ym and 'ipa respectively denote gYn and g'^a,!! ( i-^- the diagonal fluctuation of and 
9a respectively), ( ) expresses the the expectation value, and the operators O^, Om and 
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the kernel T„ ^(t, t') are given as shown below: 

O, = -iTi jdr CS,{-drA + i [S™,X„]) , (4.33) 

N+l 

Om = (C'/l^m,l/ — Cl/^m,/l) , (4.34) 

1=2 

T„,„(r, r') = SmnSir - r') - i,b^Vn{r')Orn{r)) . (4.35) 

Note that these quantum-corrected transformation laws are much more involved than 
those at the tree level. Even the linear law 6eXm = —ie^'jmQ for the bosonic field gets 
modified non-trivially contrary to naive expectation. It is not simply the expectation 
value of the linear law itself because what is relevant is the effective law that acts on F. 

4.2 Calculation of the effective SUSY transformation laws 

Actual evaluation of Ar^ and A6a at 1-loop to the relevant order is essentially straight- 
forward but extremely cumbersome. Since similar calculations performed at order 2 were 



fully displayed in the previous work[|l4|, we shall only give a sketch of the general procedure 
for the present case. 

When the variations 5^ and 6b are performed and the definitions of the composite 
operators and Om are substituted into ([4.31|) and (|4.32|) , it is not difficult to see that, 
at 1-loop, terms contributing to Ar^ and A6a can be expressed as Feynman diagrams 
composed of products of tree-level 2-point functions for various fields. These are of course 



computed from the quadratic part of the gauge-fixed action already displayed in ( |3.16|) 



The only complication is that we must disentangle the mixings among fields, which, due to 
the presence of the fermionic background, includes those between bosonic and fermionic 
fields. The results for the non-trivial 2-point functions are given in the Appendix B. 

Although most of these diagrams produce local expressions as desired, but there exist 
some which give non-local contributions. For example, a term contributing to Ar^ is of 
the form 

/ [dr'dr" {Cir)C\T")){Cuir")C^,ir')){(-A + ^r,^Y^ (/')>;; (r))e^^^(r') , (4.36) 



where the second factor (Cii(r")Cii(r')) going like 1/d produces unwanted non-locality. 
Fortunately, when one isolates all the terms of this sort, one notices that they precisely 
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cancel in the SUSY Ward identity due to the following BRST Ward identity: 

{5BQ{r")) = - jdr (^^;:^(^By.(r)Q(r")) + ^(5BV^a(r)Q(r")) j , (4.37) 
Q^r") = -iC{T") (-1* - ir„y;) (r") + iC*{T") (-1 + ir„y„) (r") . (4.38) 

We omit the details of the demonstration. 

Just as for the effective action, the calculation of the relevant diagrams for the SUSY 
transformation laws was performed with the aid of Mathcmatica. Even after various 
simplification procedures the results are quite complicated, with the longest expression, 
A9a with 4 ^^'s, consisting of 56 terms. So we relegate them to the Appendix C in order 
not to interrupt the flow of the main text. A glance at this result would convince one 
that, when we prove the invariance of efi^cctivc action under these transformations in the 
sequel, calculations of both the effective action and the SUSY transformations must be 
correct. 

4.3 SUSY invariance of the effective action 

Now we come to the last stage of this work, the check of the invariance of the effective 
action under the SUSY transformation laws just computed. Logically what we wish to 
demonstrate is quite simple: The tree level SUSY variation of the 1-loop effective action 
and the 1-loop level SUSY variation of the tree action must cancel. The difficulty is 
again that above 0{9'^) we need judicious applications of various Fierz identities to effect 
such cancellations. The required Fierz identities are more involved than those used in 
the process of simplifying the effective action itself, since we have one more independent 
spinor, namely the SUSY transformation parameter e^. Below we give a sketch of how 
such cancellations take place, choosing the case involving r^4^ as an example. 

For the structure with one e and 3 ^'s, what we wish to prove is 

(5r),,,3 = 5f )f + +si'^f^;} + <5Wr(°) + 5(^)r(°) = o . (4.39) 

Substituting the SUSY transformation laws and bring the result into certain standard 
forms by integration by parts, ST turned out to contain 40 terms. They arc classified into 
two groups, one with 3 free indices and the rest with 1 free index. We then apply three 
independent 3-frec-index Fierz identities to reduce the terms in the first group down to 
1-free-index type. Since the identities used are somewhat complicated, let us display only 
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one of them as an example: 

- 2{efe) {e^^^e) + 2{ee) (ef^'^e) - 2{e^''^e) {e^^'^'^^e) 

- 2{ee) {ejW) + {er'O) {9r''9) Sik + 2{e9) Sik 

- 2{er'9) {9r''9) Sik + 2{99) (ef^) 5jk - (e^'O) {9r''9) 5jk 

- 2{e9) {9f9) Sjk + 2{eY'9) {9Y''9) 5ju ■ (4.40) 

After this reduction, 5V consists of 29 terms, all with one free index only. To this we 
apply the following relatively simple 1-free-index Fierz identities: 

(i) (67^1^) (^7^1^^) = 2{99) {eY9) - {€^''9) {9Y'9) 

- 2{e9) (9^9) + {eY^9) {9Y^'9) , 

(ii) (e7"i^') (^7"i^^) = 2{99) {ef9) - {€^"^'9) (9^'''9) 

- 2{e9) {9Y9) + (€-^9) (h^'^'O), 

(iii) {e-f^'9) {9Y''9) = -2{99) {€^9) + {e-f''''9) {9-/^' 9) 

+ 2{e9) {9f9) + (e7"i^) {9-i''^'9), 

(iv) (e7"i"^^) (^7"i"2^^) = -2{99) (ef^) - 2(e7"i^^) (^7"^^) 

- 6(e^) {9^'9) - 2(e7"i^) (^7"^^^) , 

(v) (e7"i"2^) (^7"i"2^^) = -2(e7"i^) (^7«i'^) . 

Then (5r)e_g3 can be brought to the following form, in which the 4 fermions involved are 
all different: 

= y^^ \ — 167^ — + — 167^ 167^ — 

bivj, (e7"^"^6'') (^7"i"^'^g) bin, {e9) {9Y'9) bin, {e9) {9f^9) 
5iri, (e7"i^) (6'7"i^i^) 5 i rj, (e7"i^) (^7"i^i^) \ 

^ 167 ^ 167 ; ■ ^ ^ ^ 

Finally, by a Fierz identity of the type described in the Appendix A, this vanishes iden- 
tically. 
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The vanishing of {Sr)^^0 can be proved similarly with much less effort. On the other 
hand, for {ST)^ g5 and {6T)^ g7 we encounter a considerable difficulty: We needed up to 5- 
free-index identities which are in general quite complicated^ and moreover for expressions 
with more than six spinors there are many possibilities to choose the four among them to 
which to apply the Fierz. If the choice of which identity to use or which spinors to act on 
is not appropriate, the result immediately becomes more complicated, by dozens of terms, 
than the one prior to the application. For these reasons, even after considerable amount 
of trial end error efforts we could not fully identify all the Fierz identities responsible for 
the vanishing of the variation. To overcome this difficulty, we took the following strategy: 
First simplify the expression as much as possible by the judicious application of various 
Fierz identities. Then, when it gets reduced to a sufficiently simple form, check if it 
vanishes "numerically" by the use of explicit representations of the 5*0 (9) 7-matrices. In 
this way, we succeeded in checking the full SUSY invariance of our effective action. 



5 Summary and Discussions 

In this paper, as a necessary step toward clarification of the power of supersymmetry in 
Matrix theory, we have performed a complete calculation of the 1-loop effective action 
for arbitrary off-shell trajectory and spin degrees of freedom for a probe D-particle at 
order 4 in the derivative expansion. Although the calculation was quite involved in the 
intermediate stages, the result presented in ( |3.26|) ~ ( p.29| ) turned out to be remarkably 



simple. Further, for an independent check of this result as well as for its own interest, we 
have computed the quantum-corrected SUSY transformation laws to the appropriate order 
based on the SUSY Ward identity developed previously |]14|- This computation was again 
extremely cumbersome and we obtained rather complicated results shown in Appendix C. 
To test the invariance of the effective action under these transformations, it was crucial 
to apply a series of judicious Fierz identities, many of which had not been known. We 
developed a new efficient algorithm and generated the necessary identities. With the aid 
of these identities and some numerical computations, we succeeded in checking the desired 
invariance. 

An obvious important question about our result for the effective action is whether 
it agrees with the supergravity calculation. Unfortunately, at present time this question 
cannot be answered for the following reason. A relevant supergravity calculation was per- 
formed in ||20| up to 0(6'^), adapting a technique developed for the case of 11-dimensional 



supermembrane in [R3] to the same order. Subsequently it was argued in [19| that by 



^The longest 5-free-index identity we used consists of 109 terms. 
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appropriate re-definitions of tlie fields r^iT) and 9a{T) such an effective action can be 
made to agree with the Matrix theory calculation to that orderQ In order to check our 
complete result including all the spin effects at order 4, one needs a calculation on the 
supergravity side up to 0{6^). Judging from the existing calculation at (9(6'^) level|32 



which contains some ambiguities left unresolved, this appears to be quite a difficult task. 
Thus, assuming that the agreement with the supergravity calculation would persist, our 
result stands as a prediction (up to field re-definitions) until such a challenge will have 
been met. 

We conclude by recalling the prime motivation that prompted this work, namely the 
study of the power of supersymmetry in Matrix theory. What is most curious is to see 
how much of the highly non-trivial yet remarkably simple structure of the completely off- 
shell order 4 effective action we computed is determined by supersymmetry alone. This 
requires enumerating the most general form of the effective action at this order and study 
how the coefficients of independent structures are restricted by the requirement of SUSY 
invariance. Such a work is now underway p3| and we hope to communicate the result 
elsewhere. 
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Appendix A: A new efficient algorithm for generating SO (9) 
Fierz identities 



In this appendix, we describe a new efficient algorithm for generating the 5*0(9) Fierz 
identities and make a remark on the existence of a class of identities which are often 
overlooked. 

Let us define two types of four fermion structures /™ and g^^, with n indices contracted 
and m indices free, as follows^ : 

/;;" = (A[7'^i-«"^i-i*A2)(Aj7'^"-''^*'=+i-*™A4) , (A.l) 
= (A^^'^i-'inn-i*A3)(Af7'^"-"i*'=+i-*™Ai) . (A.2) 

The generic Fierz identities relate these two types of fermion bilinears. (In general, some 
of the free indices may be carried by factors of Kronecker (5's such as 5*^*^, etc.) In |]^ a 
systematic procedure to generate such identities was described. A part of this algorithm 
requires generation of tensor product identities for 7-matrices and this turned out to be 
progressively time-consuming as the number of free indices increases. Besides, the set 
of Fierz identities so obtained are highly redundant, including the repetition of relations 
already obtained for lower number of free indices. 

The basic idea of our new algorithm is to make full use of the Fierz identities for m-free- 
index in the calculation of the (m + l)-free-index case, i.e. it is an inductive algorithm. 
Suppose we already have all the identities with up to m free indices. Concentrating on 
the 7-matrix structure, we have m + 1 /-type structures of the form 

^ai...a„n...ifc^an...aiifc+l...jm 

expressed in terms of the g-tjpe structures. Now we try to add another free index j to 
this relation. Clearly there are 4 places to insert 7-'. So we generate 4(m + 1) relations. 
This is considerably smaller than 4'""'"^ relations generated by the method of [0. 

For example, adding 7-^ to the left-most position, we get, in an obvious tensor product 
notation, 

^J^ai...a„ji...ife ^a„...aiik+i...ini 

= ^_]_^™^ai...a„jii...ife ^ ^a„...aiife+i...im 
_j_ ^^a.i...a„-iii...ik ^ ...ai 

_|_ ^ ~ ' ^ ^_^yi+l-l^jii^ai...anii..-ii---ik (g, ^an---ai«fe+i---«m ^J^ 



'Labels which distinguish the various ways the free indices are distributed are suppressed. 
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where ii means that it is deleted. The first term contains m + 1 free indices, while the 
second sum contains structures with m — 1 indices, which can be reduced into g-tjpe 
structures by the relations already computed. 

Now on the right hand side of the original m-free-index Fierz relations, we have g-tjpe 
structures of the form (** stands for a set of indices) 

J2^irclTa- (A.5) 

When is added in the manner above, this becomes 

J2circh**i'u, (A.6) 

which can easily be computed. This produces g-tjpe structures with one more or one 
less free indices on the 7's. Equating the left and the right hand sides, we produce a 
(m + l)-free- index identity. 

Now we make a cautionary remark in enumerating all possible Fierz identities when all 
the spinors involved are distinct. The Fierz identities discussed so far (and in 0) are the 
ones where the order of the spinors are cyclically rotated. Schematically, the relation is 
of the type (Ai7A2)(A37A4) (A27A3)(A47Ai). Further cyclic rotation of course does not 
produce new identities. However, there is an additional structure (Ai7A3)(A27A4), which 
cannot be reached by cyclic rotations from the original. Therefore, one must also add 
the transformation of the type (Ai7A2)(A37A4) (Ai7A3)(A27A4). One should further 
be aware that in general these two classes of identities may contain redundancy. So the 
correct procedure is to generate all the identities of both classes and then re-solve them 
as coupled equations to find the truly independent and complete set of identities. For 
example, the identity which states the vanishing of ( [4.41| ) was obtained only through this 
procedure. 

In this way, starting from the 0-free-index identities given in we have generated all 
the independent Fierz identities up to and including 4 free indices and a part of 5-free- 
index ones. Unfortunately, the result is too space-filling to be displayed in this article. 

Appendix B: Two-point functions needed for the calculation 
of SUSY transformation laws 

We list the non-trivial two-point functions which serve as the basic elements in the cal- 
culation of SUSY transformation laws. In the expressions below, and $m are the 
(1 + 9)-component fermionic and bosonic vectors introduced in (p.l4|) , and D^\jj^ and 
bosonic and fermionic propagator matrices described in (|3.19|) and ( |3.21|) 
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respectively, with the component indices M, N etc. (running from to 9) displayed ex- 
plicitly for clarity. Further, and D^^, which were not explicitly given in the main 
text, are defined as 

t-i _f 2iv^ 

_/ A(l - 4t;^Aw^A)-i -2iA{l-4:V^Av^A)-^v"'A\ 
~ \2iAv''A{l-Av^Av^A)-'^ A(5„^ - 4t;"At;"^A)-i J' 

Dl-' ^{d- f.j)-' = -{d + - A^,,)-'A • (B.2) 
With these notations, the basic two point functions are given by 

($^(T)$^,(r')) = {r\{5ML - D^l^EK,D^'^,E,r.r'D^l^\r') , (B.3) 

(vl/;(r)$M(r')) = (r| - (^m^ - D^UEK,D^;;^E,L)-'D^lpEp,D^l\r') , (B.4) 

(*;(t)*^(t')) = {r\{5^x - ^?i5pKi^Bk5LA)-'^pi|T') , (B.5) 

($^(r)v^^(r')) = (r| - - Z^p^^S,^Z}Bi^S^,)-^Dp,V2.p^BPivk') , (B.6) 
while their "conjugates" take the form 

i^MirWAr')) = {r\{6ML - ^y;.5i.p^^;^5aL)-^^^-/^|T') , (B.7) 

(^«(r)$l,(r')) = {r\(5ML - L'y;,S^p4;^S.L)-^4LpSpA4Aik') , (B.8) 

(*.(r)*;(r')) = (r|(5,A - 4i5p/.^lKL5LA)-'^iA>') , (B.9) 

($^(r)*;(r')) = (r|(5„, - L>t;iH^^£,W^H^^)-i£,t-i5^^^J-i^|^/^ _ (B IO) 

Of course in the actual calculation, we must expand them to the appropriate order in the 
derivative expansion. 



20 



Appendix C: Quantum-corrected SUSY transformation laws 



In this appendix we display the 1-loop corrections to the supersymmetry transformation 
laws needed for checking the invariance of the order 4 effective action. The results are 
classified according to the number of ^'s in excess of the tree-level laws. 

Corrections at 0(9^): 

2r'^ Ar'^ 8r^ 

5ig^{r-a) {e-f'^O) big'^{r-v) (£7'"^) i g'^ {e-f"'9) 

4r^ 4r^ ' v ■ y 

35 i g^ Vi (r • vf {eY)^, ^5ig^Vi{r-a) (eY)^ 
Corrections at 0{9'^): 

^^m^_ 15 {96) (e7"'g) 3big^r^,n, je^e) {OfW) _ lObi g^ ie^e) jeyW) 
32 64 r9 64 

35 i g^ Vi Vj {e-f"'9) {Of ^6) 35 ig'^{r- v) {efO) {Of "^9) 5 i g^ {e^Q) {Q-f'^Q) 
^ 6479 ^ ^4^9 ^^77 

_ l^^ig^r^Vj je-f^e) {Of "^9) SBig^nvj je-f^O) {O^^'^e) _ 105 ig^nvj {eO) {Of ^"^6) 
64r9 ^ 64r9 64r9 

15 ig' {e9) {O^^e) _ 5ig\e^'d) {d-f'^d) 

32 32 ■ ^ ' ' 
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175 1 n Vj Vk {OfH) (7^e)^ big-'v, {Of ^9) (7^e), 35 z n (r ■ v) {ef'O) {l'^)a 
64 r9 8r7 ^ 16 

32 16 32 

15 i g^ Vi {9f^^9) (7^'^e)^ 35 i g^ u (r • v) {9f^^9) (7^'^e)^ hig'^n {9Y^''9) (7^'^e)^ 
^ 647^ 3279 ^ 647^ 

Big^Ti {9Y^''9) (7^'^e)^ 105 i g^ u Vj vu {9Y''^9) (7^'^e)^ 105 i g^ u {ef9) 9^ 
^ 327^^ ^ 6479 ^ §^79 

105 i g^ tti {ef9) 9^ ^ 245 i g^ Vi (r ■ v) {ef9) 9^ 315 i g^ u (r ■ vf {ef9) 9^ 



32 16 r9 8r 
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35iffV,(r-a)(eyg)ga 15iff^t;, (eY^) 35 i g V, (r ■ tQ (ef g) 

5iy^r,(ey^>^ 105iff^^V,(eg)(yg), 25 i ff^ a, (e^) (fg), 
16 16 r9 32 

_ 175 i g" Vi (r • v) {e9) {f9)^ ^15ig^ {e9) {f9)^ _35ig^ r, (r • v) {e9) {f9)^ 



+ 



32 r9 32 16 

5 7 (r/9) 105 7 (/^ r, (n ''^) (7'^^), 5 , </ a, {n "0) iy'O)^ 



16 64 r9 32 

,2 .. f .. .:\ I iin\ I . in\ 1 nr. : .:2 ..2 



32 r9 64 r9 16 

35 i g^ (r • v) (e^^9) (Y^9)^ big'^Vj {ej^9) {Y^9)^ 35 i r, (r • t;) (e7^g) (7^^^)^ 

16 r9 32 16 

5ig^n je^'O) {f'0)a 105 ^ n Vj vu {ti^9) {-^H)^ hig'v, (67'^) 

16 32 4r7 



16 r9 16 4r7 

35 % g" n (r • v) {e9) {f9)^ ^5ig^u i^d) {YO)a ^^9^v^ {^I'O) {f^9) 



+ 



16 r9 16 16 

35 1 g" n {r ■ v) iey'9) {f^9)^ big'^n {ey^9) (fW)^ big'n (6^9) 9^ 
16 r9 16 16 



16r7 16r7 ' ^ ' ^ 
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Corrections at 0{6'^): 

~ 128 ^ 12879 

128 128 r^i 



63 i n Tj Vk (Ol^'O) (Ol^'O) (Ye)^ lig^n (66) (efW) (7^e), 
256 128 r9 

_7ig^Vi {6Y''6) (6-f^''6) (7^e)„ 63 i g^ n (r • v) {6Y''6) (6-f^''6) (7^e)^ 
256 r9 ^ 128 

, 7>gh; (0^/-''''0) (0^,'-0) (7'r)„ 21igh; (Oy'^O) (O^/'O) (-^e)^ 



(C.5) 



+ 



+ 



32 r9 128 r9 

7ig^n {Ol'% {01'% (7^'g)a 189 i n Tj Vk {6f'6) {6j^'6) (j^e)^ 
64 r9 256 

315 i g^ n rj {6^% {6-f^% (Ve)^ 189 i g" n Tj Vk {6Y'"'6) {6Y''"'6) (V e)^ 
256 256 

7ig^ri {6-f^''6) {6-i^6) {fh)^ 63 i g'^ n r-j Vk (^7'="^) (^7^'™^) (7^'e)^ 
128 r9 256 

7ig'^ri {66) {6f^^6) {i^^e)^ 21 i Vi {6-f^^6) {6f^^6) {i^'^e)^ 
32^9 256 r9 

63 i g^ ri (r • v) {6y^^6) {6Y^^6) {i^^e)^ 7ig^ri {6y^^6) {6^6) (7^'^e)„ 
^ 128 ^ 128 r9 

35 i g^ Vi {6f^d) {6y^6) {i^^e)^ 7ig^ri {6f^^6) {6y''^6) {i^^e)^ 
^ 12879 ^ 6479 

21 i n {6-f^^6) {6i'^^6) {i^^e)^ 189 i g^ n rj Vk {6Y"'6) (^7^'^™^) {^''^e)^ 
12879 ^ 256 

63 i g^ Ti Tj Vk {6Y^6) {6-f^''"'6) (V'^e)^ 63 i g^ u rj Vk {6Y''6) {6-f^^"'6) (7'"'e)^ 
64rii ^ 256 r^i 

35zgV, {66) {ef6) 6^ 7ig^v, {e-fW) {6Y'6) 6^ 
64 r9 128 r9 

63 i g"^ Ti (r • v) (e-fW) {6-f'W) _ 35 i g^ u {e6) (6^6) 6^ 
128 r^i 64^9 



+ 
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64 r9 64 r9 

189 i n rj Vk {e9) {O^^H) {^9)^ ^ 63 i ^. ^ . (^y ^) (^yfei^) (y^)^ 



+ 



+ 



64rii 32rii 

35 i g" n (e^'9) {9^^ 9) (^9)^ 2\ig^v, (e9) (9f^9) j^W)^ 
64 r9 128 r9 

63 ig'^n (r • v) {e9) {9Y^9) (^^9)^ lig'r, {e9) {9Y^9) j^W)^ 
32rii 64 r9 

lig'^u {ei'9) {9Y0) {i^9)^ lig^n {ef9) {9^^9) {^^9)^ 



+ 



64 r9 64 r9 

35 1 g" r, {e9) {9fW) {^^9)^ 35 % g" v, {e^^d) {df^'^d) {^%^ 
64 r9 128 r9 

63^5^^ (r • v) {e-f^9) {9^'^^9) (7^^)^ _7ig^ri {e^^9) {9^'^^9) {^^9)^ 
32rii 6479 

2,hig^ri{e^W){9f^%{^%^ 1 1 g^ n {99) {e^^ 9) {f^ 9) ^ 
64 r9 16 r9 

lig^n {e9) {9-i^9) {f^9)^ lig'v, (67^'^) (^T^'^'g) (y^^)^ 
16 r9 128 r9 

lig'ri{n'9){9^^%{:fH)^ 6'iig'r,r^v^{e^'9){9^H){:fH)^ 
64 r9 32rii 

63 i g^ n rj Vk {e-f^9) {9-i^^9) {Y^9)^ 189 i g^ n rj Vk {e-i^9) {9^^^9) {f^9) 



+ 



+ 



+ 



128 128 

189 i g^ u Tj Vk {e9) {9-f^^^9) {Y^9)^ 63 i g^ u Vj Vk (eYd) (^7*^'^) h^^9)^ 
128 128 

63 i g'^ u Vj Vk ie9) {9-f^^9) {f^^9)^ _ 63 i n rj Vk {e-f^9) {9-f^^"'9) (7'^"^^)^ 
128 r^i 128 r^i 

35iyV, {e^^9) {9f^9) 9^ lig^n (e^) {9^'' 9) (7^^)^ 

lig'r, (e7^g) {9-^'^% {-jH)^ lig'^n {ei'9) {9-f^% {^"9)^ 
16 r9 64 r9 
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Corrections at 0(6^): 



+ 



+ 



1024 1024 

231 i(/ r, Vj r,, (O-'^'O) {e-'"""0) (-'c)„ 21 1 (/ r, {ey>'e) {0-'"'0) {e- """ff) {-^^e)^ 



+ 



1024 1024 

21 i n (e-f^^e) (e-f'^'^e) [ef^'^e) (-r^'^e)^ _ 231 i n (ef^e) (e^^'^e) (^7'="*"^) (7'"*e)^ 

2048 r^i 1024713 

231 i n rj Tk (^7*'"^) {O^^^'W) {O-i^'^W) (V^'e)^ 21 i g'^ n (eO) (e^'^e) {d-i^^O) {^W)^ 
2048 ^ 512 

23HgV,r, rfc(6g)(g7^'g)(g7^'g)(7^g), 2HgV, (e7^g) (^7^'^) (^7^^'^) (7^g)^ 
512 512 

21 i g^ n (e-f^e) (e-f^^e) (ef'^^e) (7*^^)^ 231 i / (eye) (ey^^e) (^7^'"*^) (7^=^)^ 

^ 512711 ^ 512713 

21 i g^ n (e-f^e) (e-f^^e) {e-i^^e) (f'^e)^ ^ 231 i g^ rj Tk (ee) (ef'^e) (b^^^^'B) {-i^^b)^ 



512 rii 512 ri3 

231 i g^ n Tj Tk (e-f^B) {Bf^B) {B-f^'^B) {-f^'^B)^ ^ 231 i g^ rj n (eYB) {Bf^^'B) {B-f^'^'^B) {-f^'^B)^ 



512 ri3 512 ri3 

(C.7) 
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